1 



Quark-Hadron Phase- Transit ion in an Extended NJL Model 
with Scalar- Vector Interaction 

Finite Temperature and Baryon Chemical Potential Case 



Tong-Gyu Lee, 1 Yasuhiko Tsue, 2 Joao da Providencia 3 , 
Constanca Providencia 3 and Masatoshi Yamamura 4 

1 Graduate School of Integrated Arts and Sciences, Kochi University, Kochi 

780-8520, Japan 

2 Physics Division, Faculty of Science, Kochi University, Kochi 780-8520, Japan 
3 Departamento de Fisica, Universidade de Coimbra, 3004-516 Coimbra, Portugal 
4 Department of Pure and Applied Physics, Faculty of Engineering Science, 
£NJ 1 Kansai University, Suita 564-8680, Japan 

(Received July 9, 2012) 

5 ' 

\q \ The quark-hadron phase transition at finite temperature and baryon chemical potential 

is investigated in an extended NJL model with scalar-vector eight-point interaction by com- 
i r paring the pressure of symmetric nuclear matter with that of the quark matter. As a result, 

' the extended NJL phase diagram is obtained in the temperature-baryon chemical potential 

O-l plane and the effects of the scalar- vector coupling constant G% v on the chiral phase transition 

JL . are summarized. It is shown that a quarkyonic-like phase in which the chiral symmetry is re- 

stored but the elementary excitation modes are nucleonic appears just before deconfmement 
f-H 1 in this model. 

§1. Introduction 
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| To understand the strongly interacting many-particle systems governed by Quan- 

0^ ■ turn Chromodynamics (QCD) under extreme conditions at high temperature and/or 

at high density is one of the most fascinating subjects in modern theoretical physics. 
Above all, theoretical studies of the phase transition between hadronic and quark- 
gluon matters and/or the phase diagram on the temperature-chemical potential plane 
for quark-hadron many-body systems at finite temperature and density are the most 
recent interests. In these extremely hot and/or dense environment for quark-hadron 
£> ■ systems, there may exist various possible phases with rich symmetry breaking pat- 

ternfD The extremely high temperature system which is reproduced experimentally 
by the relativistic heavy ion collisions (RHIC) has been examined theoretically by 
the first principle lattice calculations. In the finite density system, however, the 
lattice QCD simulation is not straightforwardly feasible due to the so-called sign 
problem, namely, it is difficult to understand directly from QCD at finite density. 
Thus, the effective model based on QCD can be a useful tool to deal with finite 
density system. By using the various effective models, the chiral phase transition 
has been often investigated at finite temperature and density. However, it is still 
difficult to derive the definite results on the quark-hadron phase transition due to 
the quark confinement on the hadron side. 

The Nambu-Jona-Lasinio (NJL) modeP^ is one of the useful effective models 
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of QCD. This model provides wide information for hadronic systems based on the 
dynamical chiral symmetry breaking and its restoration.® The NJL model was 
originally introduced for nucleonic/hadronic degrees of freedom.®'® Some of similar 
models describe a nucleon as a dynamical quark-diquark bound-state.® However, 
the NJL model is mostly used for quark degrees of freedom in the modern form.® 
When dealing with the symmetric nuclear matter, it is necessary to reproduce the 
property of nuclear saturation as the Walecka mode® has succeeded in describing 
phenomenologically the saturation property of symmetric nuclear matter without 
chiral symmetry in which the nucleon is treated as a fundamental particle, not as 
a composite one. In the original NJL model with chiral symmetry, however, if the 
nucleon field is treated as not a composite but as a fundamental fermion field, it 
is unable to reproduce the nuclear matter saturation property. Here, it has been 
observed that the nuclear saturation property is well reproduced by introducing the 
scalar- vector and isoscalar-vector eight-point interaction in the orig inal NJL model,® 
the nucleon being then a fundamental fermion. Thus, it is possible to consider a 
NJL-type model in which the nuclear saturation property is satisfied, as one of the 
possible models of nuclear matter .®'Q® 

In this paper, with the aim of extending our previous work^ to finite temper- 
ature and chemical potential case, we investigate the quark-hadron phase transition 
and draw the phase diagram between the symmetric nuclear matter and free quark 
phase without diquark correlation. Then, for nuclear matter, we adopt the extended 
NJL model with the scalar-vector eight-point interaction and treat the nucleon field 
as a fundamental fermion field with N c = 1 in which N c is the number of colors. As 
for quark matter, we adopt the extended NJL model with N c = 3 and treat the quark 
field as a fundamental fermion field. In the zero-temperature results in Ref. [TP]) , it 
is found that a first-order quark-hadron phase transition is obtained at finite density 
and the quark-hadron phase transition occurs after the chiral symmetry restoration 
in nuclear matter. 

This paper is organized as follows: In the next section, we briefly recapitulate the 
extended NJL model at finite temperature and baryon chemical potential for nuclear 
and quark matters following Ref. fT0|) . In £j3l the numerical results are given and the 
quark-hadron phase transition is described in this model. In 21 the extended NJL 
phase diagram with scalar-vector eight-point interaction at finite temperature and 
baryon chemical potential is presented. Also, the dependence of the scalar-vector 
coupling constant on the phase diagram is shown. The last section is devoted to a 
summary and concluding remarks. 

§2. Brief recapitulation of the extended NJL model for nuclear and 
quark matters at finite temperature and density 

In this section, following Ref llOp . the same NJL-type models including the scalar- 
vector eight-point interaction are given for nuclear and quark matters at finite tem- 
perature and density, the model parameters and the number of colors being different. 
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2.1. Lagrangian density, gap equation and pressure 

Let us start with the following Lagrangian density for nuclear and quark matters: 

where the subscript/superscript i denotes nothing but an index which represents 
the case of nuclear matter (i = N) or quark matter (i = q). Here, tpi represents 
fermion field, that is, ipN is nucleon field and ip q is quark field, and r are the Pauli 
matrices in isospin space. The first two terms are the original NJL model Lagrangian 
density. The third term is a vector- vector repulsive term with G % v . The last term 
is a scalar- vector and isoscalar- vector coupling term with G l sv . Parameters G l v and 
G l sv represent a coupling constant of four-point vector-vector interaction and that of 
eight-point scalar-vector interaction, respectively. As is well known, the pure NJL 
interaction alone is not enough to reproduce the property of nuclear saturation in 
the original NJL model. Thus, we introduce the last two terms in Eq. (|2-ip so as to 
reproduce the nuclear matter saturation properties.^ Then, in this paper, this model 
is called an extended NJL model. This model is nonrenormalizable, so that we adopt 
a three-momentum cutoff scheme in which the cutoff parameter A{ is introduced. 

Under the mean field approximation^], the Lagrangian density Cf IF and the 
Hamiltonian density Ji\ [F are obtained as 

£f F = - - J2i$a iPi + Q , 

yMF = _^. 7 . + rmipiipi + Jiiiprfipi - Ci , 

d = -G\U^)) 2 + G* # i7 Vi)) 2 + SG^^F^TVi)) 2 , (2-2) 

where 

m, = -2 [G\ - GiMil ^)) 2 ] {®M , (2-3) 
JH = 2[Gl + GlA^f] #i7>i)) (2-4) 

for nuclear matter (i = N) and quark matter (i = q), respectively. Here, the symbol 
((•••)) denotes the finite-temperature expectation value which represents thermal 
average. 

In addition, we introduce the chemical potential fi{ to deal with a finite density 
system: 

H' t = Uf F - fn4i>i 

= -iiptf ■ Vipi + rmipiipi - uliprfipi - Ci , (2-5) 

*' We make a replacement such as i/jFtp —¥ ((i/jFtp)) + (ipTi/i — ((i/jFtp))) where t^rip is bi-linear 
quantities in the fermion fields and _T are matrices in Dirac, flavor and color space. The symbol 
{(■ ■ ■}) denotes the expectation value at finite temperature or thermal average. Here, the fluctuation, 
ipr-tp — {(tpr^)}), is linearized. Namely, we consider two non- vanishing terms: ((V'jV'i)) 7^ an d 
((V';7 < Vj)) 7^ 0- (pi = ((ipii>i)) = {{^irf^i)) where pi represents the fermion number density.) 
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where /x£ is the effective chemical potential: 

^ = fM-JU = IM-2 [Gt + GlA^Af] «&7Vi» • (2-6) 
Here, the expectation values at finite temperature is given as 

(®M = n [ ,iX K - nl) , (2-7) 



with 



i7°^»=^ / + ( 2 - 



z/i = 2N)Ni , (2-9) 



(2-10) 



where i/j is the degeneracy factor in which iV* and Nj represent the numbers of 
color and flavor, and n!j_ are the fermion number distribution functions with (3 = 
1/T in which T is temperature. Here, we have eliminated the contribution of the 
occupied negative energy states from the nucleon and/or quark number density itself 
in Eq. (|2-8|) . As a result, we obtain a self-consistent set of equations, Eqs. f|2-3[) and 
(|2-6|) - (|2-8|) with (|2-1U|) . A self-consistent equation for mi in Eq. (|2-3|) is noting but 
the so-called gap equation in the BCS theory. 

The thermodynamic potential density u>i for nuclear and quark matters are de- 
fined as 

= mD) - hW)) - ±((Si)) , (2-n) 



where 



and 



mt IF )) = ((Ad ■ p)A)) -^((AAf 

+G;((^7 ^)) 2 + Gi v ((AA)) 2 ((Al°A)) 2 , (2-12) 

«M» = ((Al°A)) , (2-13) 

((Si)) = - v *Jjffia[< lnn V + (1 " O ln(l " n\) 

+n!_lnn!_ + (1 -n'l)hi(l - n!_) 1 , (2-14) 



■M ■ p)< v» /',- / ^ « - n!_) • (2-15) 

(2^) 3 . I p 2 + m 2 



Incidentally, we obtain the gap equation in Eq. (]2 3j) and the fermion number dis- 
tribution functions in Eq. (|2-1U|) again by minimizing u>i with respect to rrtj and n±. 
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From this thermodynamic potential density, the pressure of nuclear matter (i = N) 
and that of quark matter (i = q) are given as 



Pi 



nDKPi 



m* IF )){ P i = 0)] + miNi)) + -{{S t )) , (2-16) 



where we subtracted the zero-density expectation value of the mean field Hamiltonian 
in the vacuum. This subtraction method is the same as that used in the expression 
for energy density which is given later. We will discuss the determination of the 
realized phase by comparing the pressure of nuclear matter and of quark matter. 

2.2. Model-parameters 

For nuclear matter, the extended NJL model has four parameters: Gf , Gf , Gf v 
and Aff. These parameters are determined by four conditions at zero tempera- 
ture: the nucleoli mass in vacuum, ?ri/v(pjv=0) = 939 MeV, the normal nuclear 
density, p® N = 0.17/fm 3 , the nucleon mass at normal nuclear density, mjv(p' 

■15 



N ) 

0.6m.7v(p7v = 0), and the saturation properties of nuclear matter, Wn(p n ) = - 
MeV. At zero temperature, the fermion number distribution function in Eq. (|2T0[ 

reduces to the Heaviside step function, nf = 6>(/x[ — p 2 + m? ) , and nj = 1. Thus 
in the nuclear matter case, the gap equation at zero temperature is expressed as 



m N 



-2G 



N 



N 



1 



^sv 2 



2vr 2 



(V'ivY'iv) 



d\ P \ 



P 



p 2 + m 2 N 



(2-17) 
(2-18) 



where 



PN 



Pf 



Qtt 2 



Pf 



m 



N ' 



(2-19) 
(2-20) 



and v N = 2NfNf with Nf = 2 and Nf = 1. The symbol (• • • ) denotes the 
expectation value at zero temperature. The energy density per single nucleon at 
finite baryon density and zero temperature is evaluated as 



W N {p 



N) 



PN 



m N (p N = 0) 



(2-21) 



where 



{-Hn F ){pn) = ^ n (i-p)^n) 



G 



N 



d\p\ 



^sv 2 



\P\ 



p 2 + m 2 N 



(^Af + G^pj, , (2-22) 



N 2 



(2-23) 
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Table I. The parameter set for nuclear matter (i — N). 



A N [MeV] 


1 r N a 2 1 r N a 2 1 


, ~'sv /i -N 


377.8 


] 19.2596 | 17.9824 | 


-1069.89 


Table II. 


The parameter set for quark matter (i = q). 


A q [UeV] 


| G q A 2 | qia 2 


1 G q A s 


653.961 


] 2.13922 | 


| free 



Here, we take the form of subtracting vacuum value for (^HM ) in Eq. (|2-21[) . The 
values of model-parameters for nuclear matter are summarized in Table I. 

As may be seen in this table, the momentum cutoff An is rather small. It is 
possible, however, that the cutoff similarly increases with density by considering 
a chemical potential dependent cutoff. 11 ' Under these model-parameters, a rather 
reasonable value of the incompressibility is obtained numerically as K ~ 260 MeV in 
which the incompressibility of nuclear matter at normal nuclear density is evaluated 
as K = 9p^ r d 2 M / jv(pw)/^PAr|p=pO • Although we fix the value m^{p° N ) in this paper, 
the incompressibility K at normal nuclear density is capable of becoming an input 
parameter instead of the nucleon mass mjv(p^r) at normal nuclear density since 
m N{p%) has influence on K. 

For the quark matter, there are three parameters, that is, Gi,Gi v and A q , in 
the extended NJL model with scalar-vector eight-point interaction. Here, we put 
G% = for the quark matter since G l v and G l sv are introduced so as to reproduce the 
nuclear matter saturation properties for the nuclear matter. However, we introduced 
the parameter G q sv . The reason why is that G q sv has influence on the chiral phase 
transition at finite density and temperature. The scalar-vector attractive interac- 
tion G% v makes the chiral condensate strength stronger. Namely, the chiral phase 
transition point is pushed to higher-density side with increasing G q sv . In this paper, 
we treat G q sv as a chiral phase transition tuning parameter. The parameters G q s and 
A q are determined by two conditions: the vacuum value for the dynamical quark 
mass, m q = 313 MeV, and the pion decay constant, f w = 93 MeV. For the value 
of Gi v , there is no criterion to determine it in this stage. Thus, we treat G q v as 
a free parameter. The physical quantities of quark matter case are obtained from 
the corresponding ones for the nuclear matter case with i = N — >■ i = q. Namely, 
m N -> m q , Mat -> and u N = 2NfN* ^v q = 2N q N q with Nj = 2, N? = 3. The 
value of model-parameters for quark matter are summarized in Table II. 

§3. Numerical results 

3.1. Gap-solutions and chiral phase transition 

In the quark matter case, the coupling constant for the scalar- vector and isoscalar- 
vector eight-point interaction, G q sv , is a free parameter as was already mentioned 
in previous section. Hence, we put G q v A^ = —68.4 so as to realize the chiral 
phase transition in the reasonable point. Then, we take the parameter G q v as 
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Fig. 1. Left: The dynamical quark mass m q as a function of the quark chemical potential fj, q with 
G% V A% = -68.4 at T = 0, 20 and 40 MeV. Right: The quark number density divided by normal 
nuclear density p q /p% as a function of the quark chemical potential p q with Gs V A q = —68.4 
at T = 0, 20 and 40 MeV. The solid curves and branch lines are stable solutions. The dashed 
curves are unstable solutions. 

m q(Pq/3 = P%) = 0.625m g (pq = 0). The numerical results which will be pre- 
sented correspond to this parameter set and to the gap solutions in Eq. ()2-3p with 
Gl v Aq = —68.4 at finite temperature and quark chemical potential. 

The left side of Fig.[T]shows the dynamical quark mass as a functions of the quark 
chemical potential at T = 0, 20 and 40 MeV. For zero temperature, the obtained 
vacuum value for the dynamical quark mass is p q = 313 MeV. As may be seen in 
this figure, the region with the multiple solutions for the gap equation shrinks with 
increasing temperature. 

The right side of Fig.[T]shows the quark number density as a function of the quark 
chemical potential at T = 0, 20 and 40 MeV. Here, the quark number density is given 
as a multiple of normal nuclear density p a N in the vertical axis. The relation between 
the quark number density p q and the quark chemical potential p q is obtained from 
Eq. (|2-8p . The solid curves and branch lines show the stable solutions in this figure. 
Here, the solid branch lines represent the solutions of quark number density with 
massless solution, p q (m q =0). The dashed curves represent the unphysical region. 

The gap equation has multiple solutions in certain region, so that there exists an 
unphysical region with unstable solutions there. Thus, we determine the physically 
realized solution by calculating the pressure p q in Eq. (|2T6p where i = q for quark 
matter. The physically realized solution corresponds to the largest pressure at each 
temperature. Figure [2] shows the pressure of quark matter as a function of the quark 
chemical potential. 

In the case of T = MeV, the lower density solution is realized from p q = 313 
MeV to p q ~ 326 MeV in Fig. [1] on the left where the pressure is the largest value. 
Above p q ~ 326 MeV, however, the massless solution becomes physically realized. 
Namely, the chiral phase transition for T = occurs at p q ~ 326 MeV as seen in 
Fig. [2j In the quark phase, the phase with dynamical quark mass (chiral broken 
phase) is realized for the lower quark chemical potential (p q < 326 MeV) and the 
phase with massless quark (chiral symmetric phase) for the higher one (p q > 326 
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Fig. 2. The pressure of quark matter p q as a function of the quark chemical potential fi q with 
GivAq = -68.4 at T = 0,20 and 40 MeV. The solid lines represent the pressure with gap- 
solution, while the dash-dotted lines represent the pressure with massless solution. The chiral 
phase transition at T — 0, 20 and 40 MeV occur at fi q « 326, 323 and 318 MeV, respectively. 



MeV) as may be seen from Fig. [21 In the region from p q ~ 0.38/9^ (pb ~ 0.13p^ 
where ps represents the baryon number density) to p q ~ hAlp° N (pb ~ 1.80p^) in 
Fig. [Hon the right, a first-order chiral phase transition is realized and the coexistence 
of quark phases occurs. 

In the case of T = 20 MeV, the low density solution is realized up to p q ~ 323 
MeV and the massless solution becomes physically realized from p q ~ 323 MeV 
in Fig. [1] on the right. From Fig. [21 it is seen that the chiral phase transition with 
T = 20 MeV occurs at p q ~ 323 MeV. In this case, the chiral broken phase is realized 
at p q < 323 MeV and the chiral symmetric phase at p q > 323 MeV. Here, in the 
region from Pq ~ 1.30p% (p B ~ 0.43p^) to p q ~ 5Alp° N (p B ~ 1.80p^) in Fig. □ 
on the right, a first-order chiral phase transition is realized and the coexistence of 
quark phases occurs. 

In the case of T = 40 MeV, the density solution with gap-solution is realized 
in all points (from p q = to p q ~ 318 MeV) and the massless solution becomes 
physically realized from p q ~ 318 MeV in Fig. [T]on the right. In Fig. [2] with T = 40 
MeV, two branches are smoothly connected and at p q ~ 318 MeV a chiral phase 
transition occurs, or at p q ~ 5.78p^- (ps ~ 1.93p^) in Fig. Q] on the right. In this 
case, the phase with chiral symmetry breaking is realized for p q < 318 MeV and the 
phase with chiral restoration for p q > 318 MeV. Unlike in the case of T = or 20 
MeV, at T = 40 MeV, there is no jump from the lower density solution to massless 
solution since there is no unstable density solution as Fig. [T]on the right with T = 40 
MeV shows. Hence, the coexistence of quark phases doesn't occur and the order of 
the chiral phase transition is not a first-order phase transition. In case that there 
exists a jump as in Fig. [1] on the right with T = or 20 MeV, it is clear that the 
phase transition is of first order. However, because there is no jump as in Fig. [T] 
on the right with T = 40 MeV and the slope of the pressure with respect to quark 
chemical potential is the same, namely, two branches are smoothly connected as is 
seen in Fig. [21 the order of the phase transition at T = 40 MeV is a second order. 
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Incidentally, this critical density of the chiral phase transition is lower than 
that of the quark-hadron phase transition as is mentioned below. Also, the chiral 
symmetry restoration occurs systematically at the lower density side, so that it has 
no actual physical consequences and no influence on the curve of the equation of 
state in the quark phase after the quark-hadron phase transition. 

3.2. Quark-Hadron phase transition 

The main object of this paper is to investigate the quark-hadron phase transition 
in the extended NJL model with the scalar-vector eight-point interaction at finite 
temperature and density. In this subsection, we present a procedure for investigating 
the phase transition between the nuclear and quark matters and show the numerical 
results. 

For the quark-hadron phase transition, we follow the same approach to determine 
the physically realized phase as it has already been shown in previous subsection. 
Namely, we determine the realized phase by comparing the pressure of nuclear matter 
with that of quark matter at finite temperature and baryon chemical potential. For 
this purpose, the condition for chemical equilibrium is demanded as 



where /U v and 3p q indicate the chemical potential per baryon. By regarding this 
condition as the one for thermodynamic equilibrium between the hadron and the 
quark phases, we derive the corresponding condition for the pressure of hadron and 
quark phases as 



From Eq. (j2- 16fl . the pressure pn and p q can be calculated for nuclear matter and 
for quark matter, respectively. 

Figure [3] on the left shows the pressure of nuclear and of quark matter as a 
function of the baryon chemical potential which is equivalent to the nuclear chemical 
potential and the triple of the quark chemical potential at T = in the case of 
Gg V Aq = —68.4. As may be seen from this figure, there is a crossing point at a certain 
chemical potential value. From this crossing point, we can determine the coexistence 
of nuclear and quark phases. Then, about ps ~ 1236 MeV, the quark-hadron phase 
transition at T = occurs in the left side of Fig. [3j Thus, the nuclear phase (hadron 
phase) is realized for the smaller baryon chemical potential (pb = UN < 1236 MeV) 
and the quark phase for the larger chemical potential (ps = 3p q > 1236 MeV). 

On the other hand, Fig. [3] on the right shows the pressure as a function of the 
baryon number density at T = where the vertical axis is shown in logarithmic scale. 
This figure is depicted by using the relation of the density and chemical potential 
shown in Fig. [T]on the right. It is seen from this figure that the nuclear and quark 
phases coexist and a first-order quark-hadron phase transition occurs in the region 
from p N {= p B ) ~ 2Mp% to p q ~ 10.9p^ (p B ~ S.63p° N ). 

In the case of T = 20 MeV, it is seen that the quark-hadron phase transition 
occurs at ps ~ 1190 MeV as may be seen from Fig. [J] on the left. Thus, the 



p N (T) = 3p q (T) 



(3-1) 



p N (p N ,T) = p q (3p q ,T) . 



(3-2) 
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Fig. 3. The quark-hadron phase transition is shown at T = in the case of G1 v /\? q = —68.4. Left: 
The pressure of nuclear matter (i = TV) and of quark matter (i = q) as a function of the baryon 
chemical potential pb (= Pn = 3p q ) at T = MeV. Right: The pressure Pi/po as a function of 
the baryon number density pB / p% at T = MeV. Here, the vertical axis is shown in logarithmic 



scale. The pressure is divided by po where po = 
given in multiples of normal nuclear density p° N . 



l.OMeV/fm . The baryon number density is 




H B [ MeV l Pb/Pn 1 



Fig. 4. The quark-hadron phase transition is shown at T = 20 MeV with G q sv A q = —68.4. Left: The 
pressure of nuclear matter and of quark matter as a function of the baryon chemical potential 
Pb at T = 20 MeV. Right: The pressure presented on a logarithmic scale Pi/po as a function of 
the baryon number density Pb/p°n at T = 20 MeV. 



nuclear /hadron phase is realized in the region of fijsr < 1190 MeV and the quark 
phase in the region of 3/iq > 1190 MeV. From the right side of Fig. HJ it is seen 
that the coexistence of nuclear and quark phases occurs from pn(= Pb) ~ 2.49/9^ to 
p q ~ 9.99/9^ (pb ~ 3.33/9^) and a first-order chiral phase transition is realized there. 
As compared with the case of T = 0, it is seen that the crossing point (quark-hadron 
phase transition point) moves toward the lower left (small chemical potential and 
low pressure side) and the density region of the coexistence of nuclear and quark 
phases becomes smaller with increasing temperature from Fig. [U 

In the process of increasing temperature, at a certain temperature, this crossing 
point vanishes. This situation is depicted in Fig. [5] with T = 40 MeV. In this figure, 
there is no crossing point which represents the transition point, so that this leads 
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H B [MeV] p B / Pn' 



Fig. 5. The first-order quark- hadron phase transition disappears at T = 40 MeV in the extended 
NJL model with G* av A\ = -68.4. Left: The pressure pi as a function of the baryon chemical 
potential \ib at T — 40 MeV. Right: The pressure presented on a logarithmic scale Pi/po as a 
function of the baryon number density Pb/p% at T = 40 MeV. 



to the situation in which the first-order quark-hadron phase transition has already 
finished. 

§4. Phase diagram 

In this section, we present the phase diagram in the extended NJL model with 
the scalar- vector eight-point interaction at finite temperature and chemical potential. 
Then, in the quark matter, we also investigate the effects of the scalar- vector coupling 
constant G q sv which has influence on the chiral phase transition in the phase diagram. 
Namely, we discuss how the chiral phase transition line is affected by varying the 
strength of the scalar-vector interaction. 

4.1. Phase diagram with G q sv = 

First, we present the phase diagram without scalar-vector interaction. Thus, for 
the quark matter, we investigate the chiral and quark-hadron phase transition by us- 
ing the original NJL model Lagrangian density without the scalar-vector interaction, 
that is, G q sv = 0. 

Figure [6] on the top shows the phase diagram with G q sv = as a function of 
temperature and baryon chemical potential. Here, the vertical and horizontal axes 
represent temperature T and baryon chemical potential [J>b {= UN = 3/x g ), respec- 
tively. The solid curve represents the critical line of the first-order chiral phase 
transition and the dotted one represents that of the second-order chiral phase tran- 
sition. As may be seen from this figure, the critical line of the first-order chiral 
phase transition emerging from a point in the T = and /xb ~ 978 MeV terminates 
at ([ab,T) ~ (842,80) MeV. Also, the critical temperature at vanishing chemical 
potential, /jlb = 0, is found to be about 190 MeV. The dash-dotted curve represents 
the first-order quark-hadron phase transition. In this figure, it is seen that there is 
the endpoint of the first-order quark-hadron transition at (//b,T) ~ (955,39) MeV. 
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Fig. 6. The phase diagrams in the extended NJL model with G q sv —Q (top), G q sv A s q = —68.4 (mid- 
dle) and G% V A^ = —81.9 (bottom) are depicted as a temperature-baryon chemical potential 
(T-hb) plane. The solid, dotted and dash-dotted curves indicate the lst-order chiral phase 
transition, the 2nd-order chiral phase transition and the lst-order quark-hadron phase transi- 
tion, respectively. The endpoint of the lst-order chiral phase transition is at (/j,b,T) ~ (842, 80) 
MeV (top), (fi B ,T) ~ (964,26) MeV (middle) and (fi B ,T) ~ (979,1) MeV (bottom). For the 
lst-order quark-hadron phase transition, the position of the endpoint is at ((j,b,T) ~ (955, 39) 
MeV. 



4.2. Phase diagram with Gi v A^ = —68.4 

Next, we show the phase diagram with a scalar-vector interaction. For the 
scalar-vector interaction strength, we consider the value G q sv A s q = —68.4 used in the 
previous section. 
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Figure[6]in the middle shows the phase diagram in the T-pB plane with Gl v A q = 
—68.4. The critical line of the first-order chiral phase transition is depicted as a solid 
curve and that of the first-order quark-hadron phase transition as a dash-dotted 
curve. The dotted curve represents the critical line of the second-order chiral phase 
transition. This critical line of the first-order chiral phase transition emerges at T = 
and ps ~ 979 MeV and terminates at (ps,T) ~ (964,26) MeV. For the critical line 
of the first-order quark-hadron phase transition emerging from a point in the T = 
and hb ~ 1236 MeV, it terminates at (ps,T) ~ (955,39) MeV. Then, the endpoint 
of the first-order quark-hadron phase transition is located on the critical line of the 
chiral phase transition as may be seen from Fig. [6] in the middle. 

Here, it should be noted that there is a region where the quark-hadron phase 
transition occurs after the chiral symmetry restoration in the nuclear phase side. 
Namely, this suggests that a phase which is chiral symmetric but the elementary 
excitation from there is nucleonic (hadronic) could exist just before the phase tran- 
sition from the nuclear phase to the quark one. Recently, McLerran and Pisarski 
have proposed a new state of matter, the so-called quarkyonic matter,^ which is a 
phase characterized by chiral symmetry restoration and confinement based on large 
N c arguments. Thus, this chiral symmetric nuclear phase predicted by our model 
may possibly correspond to the quarkyonic phase. 

4.3. Phase diagram with Gl v A q = —81.9 

Finally, we show the phase diagram with the stronger scalar-vector interaction. 
We set the value Gi v A q = —81.9, in which m q (p q /3 = p^) = 0.63m q (p q = 0). 

Figure[6]on the bottom shows the phase diagram in the T-pB plane with Gi v A q = 
—81.9. The critical line of the first-order chiral phase transition depicted as a solid 
curve emerges at T = and ps ~ 979 MeV and terminates soon at (ps,T) ~ (979, 1) 
MeV. As for the critical line of the first-order quark-hadron phase transition depicted 
as a dash-dotted curve, the same diagram as previous two ones is obtained. 

The reason why the behavior of the first-order quark-hadron phase transition 
line seen in Fig. [6] does not change is that the critical line of the first-order quark- 
hadron phase transition exists in the chiral symmetric phase, in which m q = 
and {{^ q ip q )) = 0, that is, the quark-hadron phase transition occurs after the chiral 
phase transition for the unstable quark phase with G q sv = 0, G% v A q = —68.4 and 
Gi v A q = —81.9. Here, in the expression of the pressure p q in Eq. (|2-16p . there is no 
Gf„-dependence since {{T~i 9 MF )) in Eq. (|2-12|) and p r in Eq. (|2-6|) do not depend on 
G q sv due to U q ip q )) = 0. 

From Fig. O by varying the strength of Gi v , it is seen that the critical line of 
first-order chiral phase transition shrinks with increasing G q sv . Namely, G q v acts to 
move the endpoint of the first-order chiral phase transition toward larger baryon 
chemical potentials and the lower temperatures. In the case of Gi v A q = —68.4, the 
endpoint of the first-order quark-hadron phase transition is located on the chiral 
phase transition line. Also, the chiral symmetry restoration is shifted toward the 
larger baryon chemical potential. In the case of the stronger G$ v , the line of the 
first-order chiral phase transition disappears. 
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§5. Summary and concluding remarks 

The quark-hadron phase transition at finite temperature and baryon chemical 
potential has been investigated following Ref llOp in the extended NJL model with the 
scalar- vector eight-point interaction. In this model, as a first attempt of investigation 
of the quark-hadron phase transition, the hadron side was regarded as a symmetric 
nuclear matter and the quark side as a free quark phase with no quark-pair correla- 
tion. Here, the single nucleon in the nuclear matter and the single quark in the quark 
matter were treated as a fundamental fermions with = 1 and N q = 3, respec- 
tively. Then, in this model, the nuclear saturation property has been well reproduced 
for the nuclear matter side. On the other hand, for the quark matter side, there is 
one free parameter, G q sv . This model parameter, in this paper, was not fixed since 
there is no criterion to determine the value of G% v by using physical quantities in this 
stage. By introducing this parameter G q sv , the effective density-dependent coupling 
constant was obtained as G q {p q ) = G q {l — G q sv /G q s ■ p^). Hence, the G q sv term plays 
an important role in pushing the chiral symmetry restoration point to higher density 
side for the quark matter. Thus, the parameter G q v controls the chiral symmetry 
restoration point and/or the strength of the partial chiral symmetry restoration in 
the nuclear medium. As for the description of the quark-hadron phase transition, 
we have calculated the pressure of the nuclear matter and the quark matter, and 
determined the realized phase by comparing their pressures. As a result, a first- 
order quark-hadron phase transition is obtained at finite temperature and baryon 
chemical potential. Here, the end point of the first-order quark-hadron phase tran- 
sition is at (p B ,T) ~ (955,39) MeV with G q sv A s q = -68.4. This phase boundary is 
not changed even by varying the strength of the scalar-vector interaction because of 
G^-independence. As for the effects of the scalar- vector coupling constant G q v on 
the chiral phase transition, the critical line of the first-order chiral phase transition 
shrinks with increasing G q sv . Namely, G q v acts to move the endpoint of the first-order 
chiral phase transition toward a larger ps and a lower T. 

From the phase diagram in Fig. [6l it should be noted that there is an interesting 
phase where the quark-hadron phase transition occurs after the chiral symmetry 
restoration in the nuclear matter. This might appear as an exotic nuclear phase in 
which chiral symmetry is restored but the elementary excitation modes of matter are 
nucleonic. This phase may possibly correspond to the quarkyonic phased which is 
introduced as a chiral symmetric confined matter. 

In this paper, we have ignored the color superconducting phase.® However, this 
phase may exist at finite density system. Thus, a possible next challenging task is to 
investigate the phases of nuclear matter including the nuclear superfluidity and quark 
matter including the color superconducting state, that is, the nucleon pairing on the 
side of nuclear phase and the quark pairing on the side of quark phase. Further, it 
is widely believed that neutron star matter undergoes a phase transition to quark 
matter at high temperature and/or density. Thus, it is also interesting to investigate 
the phase transition between neutron star matter and quark matter. This leads to 
the understanding and development of the physics of neutron star. 



Quark-Hadron Phase-Transition in an Extended NJL Model 



15 



Acknowledgement 

One of the authors (T.-G. L.) would like to express his sincere thanks to Pro- 
fessor K. Iida, Dr. E. Nakano, Dr. T. Saito, Dr. K. Ishiguro and the members of 
Many-Body Theory Group of Kochi University for valuable comments and fruitful 
discussions. One of the authors (J. P.) acknowledges valuable discussions with Steven 
Moszkowski. One of the authors (Y. T.) is partially supported by the Grants-in-Aid 
of the Scientific Research (No. 23540311) from the Ministry of Education, Culture, 
Sports, Science and Technology in Japan. 

References 

1) See, for example, K. Fukushima and T. Hatsuda, Rep. Prog. Phys. 74 (2011), 014001. 

2) Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (1961), 345; ibid. 124 (1961), 246. 

3) U. Vogl and W, Weise, Prog. Part. Nucl. Phys. 27 (1991), 195. 
S. P. Klevansky, Rev. Mod. Phys. 64 (1992), 649. 

T. Hatsuda and T. Kunihiro, Phys. Rep. 247 (1994), 221. 

4) V. Koch, T. S. Biro, J. Kunz and U. Mosel, Phys. Lett. B185 (1987), 1. 

5) M. Buballa, Nucl. Phys. 611 (1996), 393. 

I. N. Mishustin, L. M. Satarova and W. Greinerb, Phys. Rep. 391 (2004), 363. 

6) W. Bentz and A. W. Thomas, Nucl. Phys. A 693, (2001), 138. 

W. Bentz, T. Horikawa, N. Ishii and A. W. Thomas, Nucl. Phys. A 720, (2003), 95. 

7) M. Buballa, Phys. Rep. 407 (2005), 205. 

8) B. D. Serot and J. D. Walecka, Int. J. Mod. Phys. E 6, (1997), 515. 

9) S. A. Moszkowski, C. Providencia, J. da Providencia and J. M. Moreira, nucl-th/0204047 

10) Y. Tsue, J. da Providencia, C. Providencia and M. Yamamura, Prog. Theor. Phys. 123 
(2010), 1013. 

11) R. Casalbuoni, R. Gatto, G. Nardulli and M. Ruggieri, Phys. Rev. D 68 (2003), 034024. 
M. Baldo, G. F. Burgio, P. Castorina, S. Plumari and D. Zappala, Phys. Rev. C 75 (2007), 
035804. 

12) L. McLerran and R. D. Pisarski, Nucl. Phys. A 796 (2007), 83. 

Y. Hidaka, L. McLerran, and R. D. Pisarski, Nucl. Phys. A 808 (2008), 117. 
L. McLerran, K. Redlich, and C. Sasaki, Nucl. Phys. A 824 (2009), 83. 

13) M. G. Alford, A. Schmitt, K. Rajagopal and T. Schafer, Rev. Mod. Phys. 80 (2008), 1445. 



